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ABSTRACT. Tobacco smoking effect on ARVs remains a topic under investigation. In this paper, 
a deterministic model is formulated by considering smoking interference in metabolism of ARVs 
and its effect on one’s adherence to drugs in order to assess stability of equilibrium states and 
determine how tobacco smoking affects antiretroviral drugs. Equilibrium states and effective re- 
production number Ref are computed and stability condition for equilibrium states established. 
Using linearization method and comparison theorem, analysis shows that disease free equilibrium 
is globally stable when Reff < 1 and it is unstable when Reff > 1. However, due to high smoking 
impairment effective reproduction number Res cannot be less than unity and the classical require- 
ment Reff < 1 for global stability of disease free equilibrium cannot be realized. Therefore tobacco 
smoking affects stability of disease free equilibrium. By applying logarithmic Lyapunov function, 
endemic equilibrium is asymptotically stable when Res > 1. The analysis shows that, as tobacco 
smoking interference with metabolism of ARVs increases, HIV infected T-cells and macrophages, 
and free virus record a corresponding increase and this shows that tobacco smoking decreases the 
efficacy of ARVs. To improve patient’s immune system and manage HIV epidemic and its therapy, 


integration of smoking cessation programs in HIV care services is recommended. 


1. INTRODUCTION 


Tobacco smoking remains a top health agenda among HIV infected patients whether they are 
under therapy or not. Contents of tobacco smoke have a devastating effect in HIV therapy and 
exacerbate pathogenesis of HIV in the in-vivo dynamics (23). The benefits of ARVs in reducing 
mortality and suffering among HIV infected individuals are negated by tobacco smoking (5; 19). 
In addition to causing immune system defective, tobacco smoking also interferes with metabolism 
of ARVs and causes drug interaction (2). Drug interaction reduces concentration of the drugs’ 
and cuts down their absorption which entails reduced efficacy. 


To gain insight on the effect of tobacco smoking when a HIV smoker is under therapy, we use 
mathematical modeling to study the behaviour of HIV within a host when ARVs are adminis- 
tered. Mathematical models play an important role in studying and understanding the dynamics 
of infectious diseases and their intervention strategies (13). In this study, the model for tobacco 
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smoking effect on antiretroviral therapy is formulated by considering HIV in vivo dynamics and 
smoking effect on T-cells and macrophages. Stability analysis of the model equilibrium states is 
performed to determine the qualitative behaviour of the dynamics. 


Stability analysis of the model steady states is one of the fundamental problem in mathemati- 
cal epidemiology (21) as it reveals important features for persistence or eradication of the disease. 
Understanding of these features which drive a particular disease helps researchers and health prac- 
titioners to design treatment and intervention strategies (28). Apart from designing treatment and 
control strategies, stability analysis also offers alert and precautions for disease outbreak. 


To perform stability analysis of equilibrium states to determine behaviour of the disease in consid- 
eration, different approaches have been proposed. However among the approaches, some are useful 
in analysis of local stability and some for global stability. Linearization method (11; 24) is useful 
in the analysis of local stability. Under this approach a system is linearized at equilibrium state 
to obtain a Jacobian matrix. From Jacobian matrix, we compute eigenvalues. If the eigenvalues 
of the Jacobian matrix are negative or have negative real parts, the equilibrium state is proved 
stable. However, the equilibrium state is proved unstable when at least one of the eigenvalues is 
positive or has positive real part. 


Sometimes it is difficult to obtain eigenvalues directly from Jacobian matrix. Whenever it is 
impossible to compute the eigenvalues from the Jacobian matrix, linearization method offers al- 
ternative methods to test the signs of eigenvalues without computing them. These methods include 
trace and determinant method, and Hurwitz criterion. Negative trace and positive determinant 
indicate that, the Jacobian matrix has negative eigenvalues. In Hurwitz criterion, we derive a 
characteristic equation and test the signs for coefficients. If coefficients do not change sign, nec- 
essary condition hold. However, sufficient condition depends on the degree of the characteristic 
equation. Linearization method is also used to analyze local stability of endemic equilibrium. 


Comparison theorem (4) and Metzler matrix (12) are used in analysis of global stability of disease 
free equilibrium. Comparison theorem proves infected classes are diminishing as uninfected classes 
are growing to attain the disease free equilibrium point. Metzler matrix forms two matrices from 
uninfected and infected classes. The method concludes global stability if the matrix from unin- 
fected classes has negative eigenvalues and the elements of the main diagonal in the matrix of the 
infected classes are negative. However, since global stability implies local stability the approaches 
may also be used to establish local stability. 


Lyapunov method can be used to establish stability of disease free and endemic equilibriums. 
The Lyapunov functions which are used in the analysis of disease free equilibrium are unique 
depending on the nature of the model under consideration. However, the analysis of endemic 
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equilibrium uses general Lyapunov functions which are reviewed as follows. Explicit Lyapunov 
function was constructed for the analysis of SEIR and SEIS epidemic models (15; 16). Logarith- 
mic Lyapunov function was used to analyze Lotka-Voltera systems (6). Later it was used in the 
analysis of endemic equilibrium for SIR, SIRS and SIS epidemic models (17). Composite qua- 
dratic Lyaponuv function was proposed and used to determine stability of endemic equilibrium 
for SIR, SIRS and SIS epidemic models (30). Later the composite-Volterra function was used in 
the analysis of endemic equilibrium for the model with relapse (30). 


As stability is a requirement for the model application in a real setting, this study analyzes 
stability of a mathematical model for tobacco smoking effect on antiretroviral therapy using lin- 
earization method, comparison theorem and Lyapunov function. This work is organized as follows: 
we describe the model in the next Section, analysis of the model is presented after the next Section, 
Numerical analysis and Conclusion mark the end. 


2. MODEL FORMULATION 


The model divides T-cells into five classes and macrophages into four classes. Antiretroviral 
drugs under consideration are reverse transcription inhibitors (RTIs) and Protease Inhibitors (PIs). 
In T-cells, X represents density of uninfected T-cells, X, density of smoking partially impaired T- 
cells, Xə density of HIV latently infected T-cells, X3 density of smoking critically impaired T-cells 
and X; density of HIV productively infected T-cells. For the case of macrophages, Y represents 
density of healthy macrophages, Yı density of smoking partially impaired macrophages, Y> density 
of HIV infected macrophages, Y3 density of smoking critically impaired macrophages and density 
of free virus is represented by V. 

A function A — — which decreases due to the presence of free virus (14; 3) is a recruitment 
rate for T-cells. The expression yX, + 7X3 represents smoking T-cells’ impairment rate with 
y < ņ being relative smoking impairment rates of X, and Xz due to the fact that, smoking 
critically impaired T-cells have high concentration of tobacco smoke poisonous and carcinogenic 
compounds. Parameters 61, $2 and T represent HIV infection rates for uninfected T-cells X, 
smoking partially impaired T-cells X, and HIV infection rate of uninfected T-cells X from HIV 
infected macrophages Y> respectively. Since reverse transcription in smoking critically impaired 
T-cells X3 is assumed to be spontaneous, in the presence of ARVs, HIV infects smoking critically 
impaired T-cells at a rate 63(1— f,¢) where e such that 0 < e < 1 is the efficacy of RTIs in blocking 





reverse transcription in T-cells and fi(@) = = i is a smoking effect in inducing metabolism 
wW 

of ARVs in smoking critically impaired T-cells and macrophages, w € (0,1) is the rate at which 

smoking induces metabolism of ARVs. When w = 0, smoking does not induce metabolism of 


ARVs and when @ = 1, smoking induces metabolism of ARVs at a highest rate. 
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Smoking partially impaired T-cells X; progress to smoking critically impaired T-cells X3 at 
arate p. HIV latently infected T-cells Xə due to the presence of smoking partially impaired 
T-cells, progress to productively infected T-cells X4 following successful reverse transcription at 


= 


a rate o(1 — fe) where f = e 7 models smoking effect in inducing metabolism of ARVs in 
smoking partially impaired T-cells and macrophages. Functions f and fı model relative smoking 
inducing effect in smoking partially and critically impaired cells as we assume smoking critically 
impaired cells experience high smoking inducing effect compared to smoking partially impaired 
cells. However, if ¢; denotes drugs’ reverse transcription blocking rate in T-cells and 7 adherence 
rate to treatment, then € = ¢,v. Parameter a is smoking induced mortality in critically impaired 
T-cells and jz; is a HIV induced mortality in HIV productively infected T-cells which produce 
infectious virions at a rate Nıuı(1 — fig). The parameter € = «VJ such that 0 < € < 1 is the 
efficacy of PIs in blocking production of infectious virions in T-cells, x is the rate at which Pls 
block production of infectious virions in T-cells and v is drugs’ adherence rate. 

Macrophages are recruited at a rate A. Expressions 64(1 — €1), 35(1 — fe.) and e(l — fie.) are 
HIV infection rates for healthy macrophages Y, smoking partially impaired macrophages Yı and 
smoking critically impaired macrophages Y3, where €; = C20 such that 0 < e < 1 is the efficacy of 
RTs in macrophages and ¢ is the rate at which RTIs block reverse transcription in macrophages. 
Expression vY + 0Y is a smoking impairment rate with relative impairment of Yı and Y; given 
by v and 0, we assume that v < 0 because smoking critically impaired macrophages have high 
concentration of tobacco smoke poisonous and carcinogenic compounds. Smoking critically im- 
paired macrophages suffer smoking induced mortality at a rate a;. HIV infected macrophages 
suffer HIV induced mortality at rate ô and produce infectious virions at a rate N2d(1 — f1&1) 
where €; = «k10 such that 0 < €; < 1 is the efficacy of PIs in blocking production of infectious 
virions in macrophages and «Kı is the rate at which PIs’ block production of infectious virions. 
Parameters uy and jt, represent natural mortalities for macrophages’ compartments and free virus 
respectively. 

The model assumes that RTIs and PIs have different efficacies in T-cells and macrophages (1). 
Metabolism of RTIs and PIs occurs at a cellular level. Tobacco smoking induces metabolism of 
RTIs and PIs uniformly. On HIV infection, smoking partially impaired T-cells join HIV latently 
infected compartment because their reverse transcription is not spontaneous and critically impaired 
T-cells join productively HIV infected T-cells due to their spontaneous reverse transcription. The 
model also assumes that unresponsiveness of T-cells and impairment of macrophages increase 
with smoking. Interaction of variables is demonstrated by Figure 2.1, state variables and model 
parameters are described in Tables 2.1 and 2.2 respectively. 
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FIGURE 2.1. Interaction of T-cells and macrophages with free virus and tobacco 
smoking in the presence of therapy. 


TABLE 2.1. Variables description 





Variable Description 





X Uninfected CD4+ T-cells 

Xı Smoking partially impaired CD4+ T-cells 
Xo HIV latently infected CD4+ T-cells 

X3 Smoking critically impaired CD4+ T-cells 
X4 HIV productively infected CD4+ T-cells 
Y Healthy macrophages 

Yı smoking partially impaired macrophages 
Yə HIV infected macrophages 

Y; Smoking critically impaired macrophages 
V Free virus 
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Table 2.2: Parameters descriptions 











Parameters Description 

A CD4+ T-cells recruitment in the absence of HIV 

k half saturation constant 

c rate at which HIV reduces newly produced CD4+ T-cells 

À recruitment rate for macrophages 

y rate of impairment from partially impaired CD4+ T-cells 

n rate of impairment from critically impaired CD4+ T-cells 

p progression rate from partially to critically impaired CD4+ T-cells 
o progression rate from HIV latent to actively infected CD4+ T-cells 
u natural mortality rate for CD4+ T-cells 

By HIV infection rate for uninfected CD4+ T-cells 

Bo HIV infection rate for partially impaired CD4+ T-cells 

B3 HIV infection rate for critically impaired T-cells 

Ba HIV infection rate for uninfected macrophages 

Bs HIV infection rate for partially impaired macrophages 

Be HIV infection rate for partially impaired macrophages 

V smoking impairment rate from partially impaired macrophages 

0 rate of impairment from critically impaired macrophages 

W HIV induced death rate for CD4+ T-cells 

[ly natural mortality rate for macrophages 

ô HIV induced death rate for macrophages 

lv natural mortality rate for free virus 

ai smoking induced death rate for impaired macrophages 

a smoking induced death rate for impaired CD4+ T-cells 

T HIV transmission rate by infected macrophages to CD4+ T-cells 

q progression rate from partially to critically impaired macrophages 
N, & No Number of virus released by a T-cell and a macrophage over the life time 
E&E RTIs efficacies in T-cells and macrophages 

E& £ PIs efficacies in T-cells and macrophages 

W Smoking induction rate in metabolism of RTIs and PIs 
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From the flow diagram 2.1, the model is governed by the following system of differential equations: 





dX cV 
(la) a pee AEA EAA 
dX 
(b) = = (yX1 + 9X3) X — BV Xi — (P+ w) Xs, 
dX: 
(1c) ax = Bi VX +7Y2X + bV X: — (o(1 — fe) +p) Xo, 
dX3 
(1d) ae T pXı — B3(1 — fıe)V X; — (a + u)X3, 
aX, 
(le) age a(1 — fe)X2 + 63(1 — fie)VX3 — (m + pw) Xa, 
dY 
(1f) ae A— Ball — 6) VY — (VY: + OY3)Y — uY, 
dY) 
Ge) GP = ON + OYS)Y — B= fa) VY = (4 + y)M, 
dY: 
(h) = = a(l- a) VY +b- fa) VY + Bo(1 — fier) V Ys — (5 + py) Yo, 
dt 
: dY: 
(li) SF =M- pl- fier) VYs — (an + Hy) Ya, 
dV 
ae Nypa(l — fr€)X4 + Nod(1 — fi€i) Yo — BiVX — bV Xı — B3(1 — fie)VX3 


(j) — Bal —€1)VY — pam fa A — e — frer)VY¥3 — eV, 


subject to initial conditions X(0) = Xo, X1(0) = X10, X2(0) = 0, X3(0) = 0, X40) = 0, 
Y (0) = Yo, Yi(0) = Y10, Y2(0) = 0, Y3(0) = 0, and V(0) = Vo. 


3. MODEL ANALYSIS 


In this section the region under which the solutions of model (1) are bounded is deduced and we 
compute equilibrium states, effective reproduction number and determine condition for equilibria 
stability. 


3.1. Boundedness of Solutions. 


To prove boundedness of solutions we consider T-cells, macrophages and free virus separately. 
If T, represents sum of T-cells in all compartments, we have: 


aT, 
- =A- — pT, — aX3 — Xa, 


H k+V 
t 


C dees 
Ape eat 
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Rearrangement gives the following equation: 
dT; 


2 — T, <A 

(2) di + ptt S 

whose general solution is given by 
A A 

(3) Te 4: (ro - =) a 
H H 


A A 
Considering the two cases when 7;(0) > — and when 7;(0) < —, we have: 
H H 


T; < ®,, where 
A 
, = max fno, =} 
u 


However, since T, represents the sum of all T-cells, it follows that: 


(4) 


(5) X, <. 


For the case of macrophages if M, represents the sum of macrophages in all compartments and we 
apply the same procedures as for T-cells, we have 





À 
(6) M < bu = maz M0, È}. 
Hy 

Since M is the sum of macrophages in all compartments, it follows that 
(7) Yo < Öm. 
We now consider equation (1j), so that 

dV 

ap Nim (1 — fi€)X4 + Nod(1 — fi€i) Yo — bV X — BoVX1 — B3(1 — fre)V X3 

—Pa(1—ea\)VY — (1 — fa) VY — bel — fires) VY — WV, 

dV 

ae Nypa(l — fig) Xa + Nod(1 — fidi)¥o — pe. 
Substitution of X4 and Y> in equation 

dV 
(8) Te < Nipa(l — fi) Xa + Nod(1 — fi€1) Yo — WV, 
yields solution which shows that free virus are also bounded. The solution is given by 
(9) V(t) < Yo, 
where 
_ Mead = FETO) | Nod = fi) M0) Nim = AHA | Nod — fi&i)A 
Woy = Max ; H : 
My My My My Hy 


The solutions of the model system are bounded in the region 


(10) Y = {(T,Y,M, V) € RP :<T<G,:05M<O,:0<V < Vo} 
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where T and M represent T-cells’ and macrophages’ compartments respectively. 
Any solution on the boundary of T converges and remains in the region. Existence, uniqueness 
and continuity of solutions of the model (1) hold in TY. 


3.2. Disease Free Equilibrium and Effective Reproduction Number Rf, ;;. 
The disease free equilibrium of the model (1) when there is no tobacco smoking and HIV is 
given by 
‘i A À 
(11) X (X, Xi, Xo, X3, Xa, Y, Vi, Yo, Y3, V) = —, 0,0,0,0, —, 0,0,0,0 
H Hy 
We use disease free equilibrium to compute effective reproduction number Resp in the next section. 


3.2.1. Effective Reproduction Number Rss. If we consider the infected classes in a model system 
(1) so that the new infections and transfer terms are defined by M; and N; respectively, the matrix 
M and N are defined by 





(12) M = 


According to (29), the effective reproduction number Res is given by 
(13) Rest = p (MN) 


which is the mazimum eigenvalue of the matrix MN~*. From the model (1), we define the new 
infections and transfer terms to be 


(yX1 + nX3)X 
ByVX +7Y2X + BoV X, 

0 

63(1 A fie) V X3 

(VY; + OY3)Y 
Ball — a) VY + 6501 — fea) VY + be — fren)VY¥3 
0 
Nim (1 — fié) X4 + Nod(1 — fi&i Yo 


and 


(p+ u)Xı 

(o(1 — fe) + u) X2 

(a + p)X3 — pXy 

(Hi + p)X4 — a(l — fe) X2 

(q + Hy)Yi 
(8 + Hy) Y2 

(a1 + Hy)Y3 — qYı 

Ha V 


(15) N; = 
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From equation (13), the effective reproduction number Rept works out to be: 








(16) Reps = max{ Reff, Reff} 
where 
Pe. os Ay | Anp 
netu)  p(p+p)(a+ p) 
Ref = a a , 
Hyld + Hy) Hyld + Hy) (or + Hy) 


The effective reproduction number Ress is given as the maximum of partial effective reproductive 
number due to tobacco smoking in T-cells Ry, and partial effective reproductive number due to 
tobacco smoking in macrophages Rrf2. Tobacco smoking dominates HIV in producing new infec- 
tions. The partial effective reproductive numbers Ref, and Refo can be written as the sum of new 
infections which are caused by smoking partially impaired cells and those caused by smoking crit- 
ically impaired cells. If the new infections which are caused by smoking partially impaired T-cells 
are defined by Repr and those caused by smoking critically impaired T-cells by Recr then 


Regi = Repr + Recr, 
(17) Anp 


(p+ n)a + u) 





A 
where Repr = eet ee Recor = 
Lup +u) 
Similarly for the case of macrophages, if the new infections which are caused by smoking par- 
tially impaired macrophages are defined by Repy and those which are caused by smoking critically 
impaired macrophages by Recm then 


Reppo = Repu + Recm, 


where 
ee) AV Avq 


—— and Recm = : 
[y(q + fy) Hylq + Hy) (Or + Hy) 





Repm = 


We use effective reproduction number Resp to determine stability of disease free equilibrium. 


3.2.2. Local stability of a Disease Free Equilibrium x°. When tobacco smoking impairs less than 
one cell (T-cells or macrophages) and HIV infects less than one cell (T-cells or macrophages) a 
case in which Refs < 1, disease free equilibrium is locally asymptotically stable. It becomes unstable 
when tobacco smoking impairs more than one cell (T-cells or macrophages) and HIV infects more 
than one cell (T-cells or macrophages) the case in which Resp > 1. To find condition for local 
stability of disease free equilibrium, we state and prove the following theorem. 


Theorem 1. : Disease free equilibrium is locally asymptotically stable when Resp < 1 and unstable 
when Ress > 1. 
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To prove the theorem, we linearize the model system (1) at disease free equilibrium to obtain a 
matrix 


A 
=i a 0 mn o o 0 TA g _¢_ A 
u u u k p 
o E ii 0 a 0 0 0 0 0 0 
A 
0 0 —wa 0 0 0 0 TA 0 bA 
u u 
0 p 0 -wo 0 0 0 0 0 0 
0 0 o(l—fe) 0 —we 0 0 0 0 0 
Oy à Or d 
9) J(x") 0 0 0 0 Gc. sige, ae a. «i _ a3 
Hy as Hy 
0 0 0 0 0 ZA ay > 0 me 0 
Hy Hy d 
0 0 0 0 0 0 0 —wa 0 2 
Hy 
0 0 0 0 0 0 q 0 -w 0 
A d 
0 0 0 0 dy 0 0 dy we. MA 
H Hy 
where 


dı = Niu (1 — fié), d2 = N2d(1 — fii), ds = B4A(1 — €1). 


Negative eigenvalues from matrix J(x°) will mean stable disease free equilibrium. Negative trace 
and positive determinant of the matrix J(x°) also means that the eigenvalues of the matrix J(x°) 
are negative and disease free equilibrium is stable. After identifying two negative eigenvalues — u 
and —uy in first and sizth columns, the matrix reduces to 


2 — w 0 = 0 0 0 0 0 
A 
0 —wa 0 0 0 eS 3H Buss 
m m 
p 0 —wo 0 0 0 0 0 
5 0 o(l— fà) 0 -ws 0 0 0 0 
ey) Ju(x) = 0 0 @ wr Pelee or A 0 
Hy Hy d 
0 0 0 0 0 -w4 0 z2 
Hy 
0 0 0 0 q 0 —we 0 
A a: 
0 0 0 dy 0 do o see iy 
H Hy 


If we denote trace and determinant of matrix J,(x°) by trJ, and detJ,, then trace and determi- 


nant are given by 


A l-e 
(21) trJ =wi(Repr — 1) + w3(Reup — 1) — A = Pal T i) 


—W4— WA — W6, 





— Py — Wo — We 
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and 
det J, = ® BA ! Ball — €1)A t1— Rew| (1 — Resp (1 — Rep pa), 
Hho My by 
(22) where Reg = (7 Ba(1 — €1) + Bippyws) CAN m (1 — AE) — fe) 


fly Ly W AW4W6 


Nð — fir€i) Pad — €1)A 
| My Hy Wa 





, B= waww wwz wwe. 


The expressions for Repr and Remp are given in (17) and (18) respectively. We find that 
(23) trJ, <0 if and only if Repr < 1 and Remp < 1, 


and 


det J, > 0 if and only if Reffi < 1, Reff2 <1 
BA | bal — e)A L1 
Hho Hyby 


24 
(24) and Reg < 





Non-negative determinant and negative trace represent a model system with locally asymptotically 


stable disease fee equilibrium when Rep < 1. However since Resp is the maximum of Reff, and 


Ref fo, stable disease free equilibrium will not be realized when Resp, > 1 and Reps, > 1. 


3.2.3. Global Stability of a Disease Free Equilibrium x°. Global stability analysis of a disease free 


equilibrium is usually done by Lyaponuv functions as applied in (27), (81), (25) and (26), and 


by comparison theorem as used in (4; 9) and (32). Since Lyapunov functions are not unique and 


present a challenge in construction, this work adopts comparison theorem in global stability anal- 


ysis for disease free equilibrium. 


Considering infected classes alone in the model (1), we write the system without uninfected com- 


partments by 
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aX +X) (2 — x) 


A 
Xi Ai (TY2 + iV) (= - x] 
X3 Xə H 
x! X : 
X’ X4 
25 * | =(M-N — À 
E cae onor (>-y) 
Y) Y» A Hy 
Y; Ys av (2 -¥) 
Vy’ V 0 
A À 
L+6V(X-—)+&(l-a)V(Y -— 
H [ly 
A À 
L= a +l a + bV Xa + B3(1 — fre)V Xs + B5(1 — fea) V Yı + BeV(1 — fier) ¥3.- 
y 
A 
Matrices M and N represent new infections and transfer terms respectively. Fort > 0, X < = 
u 
À 
and Y < — we see that 
Hy 
Xi Xı 
X5 Xo 
X X3 
X; X. 
(26) 4] <(M—N)] <7 
Yı Yı 
Y; Yə 
Y; Y; 
yv’ V 


Since the matrix M —N has negative eigenvalues, equation (26) represents a stable disease free 
equilibrium whereby (X,Y) > is a and (X 1, X2, X3, X4, V1, Yo, Y3, V) — (0,0, 0, 0,0, 0, 0, 0) 
as t + œ. This shows that a ae free equilibrium is globally asymptotically stable. We 
summarize this result in Theorem 2. 

Theorem 2. : The disease equilibrium x? is globally asymptotically stable when (X,Y) > 
E >) and (X1, Xo, X3, Xa, Ya, Yo, Ys, V) — (0,0,0, 0,0,0,0,0) for which Reff < 1 


3.2.4. Global stability of endemic equilibrium. By reverse condition as stated by (29) the endemic 
equilibrium Y* is locally stable. Lyapunov functions and LaSalle invariant principle are used to 
analyze global stability of endemic equilibrium Y*. Using a system (1), we define a logarithmic 
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Lyapunov function by 
(27) U2(T, M,V) = (7 — T*InT)+c(M — M*InM)+d(V — V* ln V), 


T, M and V represent T-cells, macrophages and free virus respectively and T*, M* and V* 
represent T-cells, macrophages and free virus at endemic equilibrium point. In full, function 


U,(T, M,V) is written as 


Us( TMEV), = bX SAX) OG = Xn) OG = KS) 
Ebi Rg = Xn 4b OG = ney EY hY 
ABV VF Va) ota (Yo = Vo in Vo) OG Ynys) 
Ebo = Vln Vv); 





Derivative of equation (28) with respect to time, gives 


























We _y (y_ X*)\ aX 14 (_ Xt) eG 4 (XB) em 
dt : Mf ae ”? Xi dt 3 X,) dt 
TE dX : DON AG 7 Y*\ dY 
(29) | oy a a Rij cde vy ai 
if VM L ag BVM 
ae ae °? Yj) d °” Yz) dt 

VV 

Hb ns Ce es 

10( va 


Substitution of the rate of change for each variable at endemic equilibrium, simplification and 


rearrangements give 

















(30) 
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where 
E TEENETE S 
R ( 7 z) (: 7 ae) — b4bs(1 — fie)V Xs (1 = ži) Pe Tx) 
sano (2) (1-2) -wra aae) (8) 
en awn (1- F) (1-2) - st - envy (1- £1) G-E) 
reer (: 7 y) (1 = ve) — bios V X3 (1 = 7) (1 a Ta) 
swan VCE Jom BY Ca 





= Yea, rS Bi 
warni- 2) (1-28) 


Following the approach in (22) and (20), F> is non-positive. Function 
dU: 
Fy <0 for X, X1, X2, X3, X4, Y, Yi, Yo, Y2, V > 0. The time derivative a <0 when 
X, X1, X2, X3, X4, Y, V1, Yo, ¥3,V > 0, for X = X*, Xı = Xi Xə = X3, X3 = X3, X4 = XG; 


dU. 
La = hep oS hg ea =Y, V=V*, a = 0. This indicates that, the largest invariant 


dU. 
set T for which — Disa singleton Y* which is endemic equilibrium. Using LaSalle Invariant 
Principle (18), T* is globally stable in the interior of T when Resp > 1. This result is summarized 
in the following theorem 


Theorem 3. : If Reps > 1, then the model system (1) has unique endemic equilibrium Y* which 
is globally asymptotically stable in the interior of Y*. 


4. NUMERICAL ANALYSIS 


In this section we illustrate the effect of tobacco smoking on antiretroviral therapy, and we show 
existence and stability of equilibrium states numerically. As smoking interferes with metabolism 
of antiretroviral drugs, stability of equilibrium states is assessed by considering infected and non 
infected classes. The simulation results on how smoking affects ARVs on the dynamics of HIV 
among T-cells and macrophages are discussed. We use parameter values in Table 4.1 to simulate 
the model using MATLAB (Version 7.1.0.246 (R14) Service Pack 3). 
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Table 4.1: Parameter values 























Parameter Value/Unit Source 

A 600 year“! Assumed 

k 12 (3) 

c 110 year™t (14) 

À 100 = (8) 

y 0.0025 year" (10) 

n 0.007 year“! Assumed 

p 0.785 year“! Assumed 

o 0.45 year"! Assumed 

u 0.135 year“! (8) 

By 0.00876 can (26) 

Bo 000012 ae (7) 

Bs 0.0016 a o (10) 

5, 0.0002 a (5) 

Bs 0.0004 a (8) 

bs 0.0005 n o (8) 
0.0016 year~* Assumed 

6 0.003 year"! Assumed 

Ly 0.775 year+ Assumed 

ly 0.0351 year+ (10) 

ô 0.25 year! (8) 

ve 50 year+ Assumed 

ay 0.03 year! Assumed 
0.102 year! Assumed 

T 0.000365 year! (14) 
0.38 year? Assumed 

N, 100 year~+ (26) 

No 100 year~+ (8) 

w ES) 0.1 & 0.5 Assumed 

G EIC 0.5 &0.6 Assumed 

k1 Éko 0.9 &0.92 Assumed 





In this section, we begin by plotting the general dynamics for model (1). 


Figure 4.1 shows the 


dynamics among HIV in T-cells and macrophages when smoking induces metabolism of ARVs 


at a rate 0.3. The dynamics show the existence of unstable disease free equilibrium between 0 
and 2 years. During this time interval, uninfected T-cells and healthy macrophages grow slightly. 
However, smoking compartments indicate a fast growth compared to HIV compartments which take 


over later. 
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FIGURE 4.1. Variation of cells’ populations in the general dynamics when a HIV 
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To show the endemic equilibrium, we plot HIV infected and smoking impaired classes and they 
are shown in Figure 4.2. Smoking impaired classes increase between 0 and 5 years and on their 
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FIGURE 4.2. Variation of HIV infected and smoking impaired classes with time. 
We have shown unstable disease free and endemic equilibriums, now we assess the effect of 
tobacco smoking on the antiretroviral drugs by looking into how HIV latently and actively infected 
T-cells, HIV infected macrophages and free virus behave when smoking interferes metabolism of 
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ARVs. 

In Figures 4.3 and 4.4, HIV actively infected T-cells, HIV infected macrophages and free virus 
increase correspondingly with increasing in interference of ARVs’ metabolism. This shows that 
smoking interference in metabolism of ARVs reduces drugs’ efficacy. For HIV latently infected 
T-cells as antiretroviral drugs are taken, they increase from 0 to 10 years. However from 10 years 
and above HIV latently infected T-cells decreases as smoking interference with metabolism of ARVs 
increases. This is the time when ARVs becomes ineffective due to tobacco smoking. 
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FIGURE 4.3. Variation of HIV actively and latently infected T-cells with respect to 


smoking inducing effect. 
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FIGURE 4.4. Variation of HIV infected macrophages and free virus with respect to 
smoking inducing effect. 

We study stability of endemic equilibrium in the dynamics of HIV in T-cells and macrophages by 
considering uninfected T-cells and macrophages, smoking partially and critically impaired classes, 
HIV infected T-cells and macrophages and free virus. Figures 4.5, 4.6, 4.7 and 4.8 show the 
trajectories for mentioned compartments. The trajectories for different initial conditions converge 
to a single point as time goes on. This indicate the existence of endemic equilibrium which we 
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already established earlier. Convergence for trajectories indicates that endemic equilibrium exists, 
and it is globally asymptotically stable whenever it exists. 
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FIGURE 4.5. Trajectories of uninfected T-cells and healthy macrophages. 
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FIGURE 4.6. Trajectories of Smoking critically impaired T-cells and macrophages. 
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FIGURE 4.7. Trajectories of HIV actively infected T-cells and macrophages. 
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FIGURE 4.8. Trajectories of free virus. 


5. CONCLUSION 


A mathematical model to determine the effect of tobacco smoking on antiretroviral therapy and 
assess stability of equilibrium states is presented and analyzed. The effective reproduction number 
Resp is computed by using next generation approach (29) and it is given as the marimum of partial 
effective reproduction numbers due to tobacco smoking in T-cells Resp, and macrophages Reffa. The 
disease free equilibrium is shown to be globally asymptotically stable when Rete < 1. However with 
the parameters that we have used in simulation, Ref p, and Rest, are not less than unit. Since Reff 
is the maximum of Rep, and Repf,, it will never be less than unity and the classical requirement 
of Resp < 1 for global stability of disease free equilibrium cannot be achieved. Therefore tobacco 
smoking affects stability of disease free equilibrium. 

Numerical simulation indicates that the unstable disease free equilibrium and stable endemic 
equilibrium exist. Convergence of trajectories in Figures 4.5, 4.6, 4.7 and 4.8 proves the existence 
of endemic equilibrium and its stability. Endemic equilibrium is asymptotically stable whenever it 
exists when Rese > 1. Simulation indicates that tobacco smoking reduces the efficacy of antiretro- 
viral therapy. This is reflected by increase of free virus, HIV infected T-cells and macrophages 
when smoking induces metabolism of ARVs. 
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